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Oxygen-blown biomass integrated gasification combined cycle (IGCC) plants are one of the most promis¬ 
ing options for clean energy generation with CO 2 abatement potential. However, the integrated nature 
of IGCC leads to difficult design problems. In this study, we present an advanced simulation environment 
for the preliminary design and retrofit of IGCC plants. We describe the modelling approach, the model 
validation strategy and the plant behaviour, as determined by sensitivity analyses. The simulation envi¬ 
ronment uses Pareto curves to examine various co-gasification and co-production case studies in terms 
of technical, economic and environmental performance. It serves as a decision support tool in the design 
stage, which can be used to explore ways to improve plant performance and to analyse the influence of 
raw materials and the unit’s operational parameters. The test and validation results are discussed. 

© 2011 Elsevier Ltd. All rights reserved. 


1. Introduction 

Expanding economies and climate change are challenging issues 
that require cleaner and more efficient power production plants. 
Integrated gasification combined cycle (IGCC) power plants are a 
promising alternative, since they allow electricity to be obtained 
efficiently from coal-based fuels, with lower costs and emissions 
per kWh than conventional pulverised coal-fired power plants. 
IGCC power plants are versatile for two reasons: a wide range of 
carbon-based inputs can be fed into the gasifier; and different prod¬ 
ucts can be obtained from the syngas that is produced, namely 
electricity, heat, chemicals, synthetic natural gas, liquid fuels or H 2 , 
depending on the specific layout of the gasification plant. 

According to NETL (2007), coal is the main feedstock for IGCC 
power plants followed by petroleum refinery by-products, natural 
gas, petcoke, and biomass. Despite this trend, renewable energy 
sources will play an important role in the future of the energy 
generation sector. Biomass is seen as an important resource in 
terms of supply prediction. It is the largest, most diverse and 
readily exploitable source of renewable energy (Ravindranath, 
Somashekar, Dasappa, & Reddy, 2004). It encompasses specific 
energy crops, residues and waste that are present in all kinds of 
societies, and it constitutes an affordable and abundant feedstock 
(Han, Mol, Lu, & Zhang, 2008) for large and small scale use in devel¬ 
oped and developing countries. Co-gasification refers to the use of 
two or more organic fuels as feedstock, one of which is of fossil 


* Corresponding author. 

E-mail addresses: luis.puigjaner@upc.edu, luis.puigjaner@upc.es (L. Puigjaner). 


origin (such as coal or fuel oil). The advantages of biomass use are 
its renewable nature and the fact that it is a sustainable source of 
energy. However, it also has some drawbacks in terms of density, 
grindability and low lower heating value (LHV), which might render 
higher C0 2 and SO x production per MW produced. Consequently, 
there is a limit to the amount of biomass that can be co-gasified in 
an IGCC plant. 

Examples of studies that have tackled the question of co¬ 
gasification are Zhao, Hao, and Xu (2006) and Koukouzas, 
Katsiadakis, Kariopouios, and Kakaras (2008), which used Aspen 
Plus® for design purposes, and had encouraging co-gasification 
results. A paper by Valero and Uson (2006) refers to the oxy-co- 
gasification of coal, petcoke and biomass in an IGCC power plant 
with oxygen-blown gasification. In addition, a study by Baliban, 
Elia, and Floudas (2010) introduced a mixture of coal, biomass 
and natural gas into the liquid gasification process. The results 
showed that this hybrid option could satisfy the entire US trans¬ 
port demand with very low environmental damage. In our previous 
study (Perez-Fortes, Bojarski, Velo, Nougues, & Puigjaner, 2009), we 
developed and evaluated a tool for IGCC power plant design that 
was based on commercial process simulation software packages. 
We used this tool to assess blends of coal, petcoke and biomass. In 
this paper we present an extended revision of our tool that can be 
applied to a wider range of raw materials. This version is based 
solely on Aspen Plus® models and involves different processing 
units. The use of just one simulation tool increases the calculation 
speed. 

For an IGCC power plant, the carbon removal technique (known 
as carbon capture and storage, CCS) must involve the use of a pre¬ 
combustion capture method, in which the high pressure that is 
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Nomenclature 

ASU 

air separation unit 

CC 

combined cycle 

CCS 

carbon capture and storage 

eg 

clean gas 

Q n 

total mass flow of inlet carbon in the feedstock 

CGE 

cold gas efficiency 

COS 

carbonyl sulfide 

Effcc 

energy efficiency for the combined cycle 

Effd 

efficiency related with the gas cleaning units 

Effglobal 

energy efficiency for the whole plant 

ER 

equivalence ratio 

GT 

gas turbine 

HHV 

higher heating value 

HRSG 

heat recovery steam generator 

HEN 

heat exchangers network 

HP, IP, LP 

’ high pressure, intermediate pressure and low pres¬ 
sure 

IGCC 

integrated gasification combined cycle 

LHV 

lower heating value 

MC 

moisture content 

MDEA 

methyldiethanol amine 

NGCC 

natural gas and combined cycle 

PSA 

pressure swing adsorption 

PSD 

particle size distribution 

rg 

raw gas 

SA 

sensitivity analysis 

ST 

steam turbine 

TOT 

turbine outlet temperature 

VL 

vapor-liquid 

WGS 

water gas shift 

WHB 

waste heat boiler 


maintained throughout the gas cleaning process is used to sepa¬ 
rate C0 2 by means of a physical absorption method (Kanniche et al., 
2010). The H 2 obtained by the gasification of solid organic matter is 
produced firstly in the gasifier itself, and secondly in a WGS reactor 
in which CO reacts with H 2 0 to produce H 2 and C0 2 . The produc¬ 
tion of H 2 inherently involves the production of C0 2 ; thus, final C0 2 
disposal solutions need to be researched. 

IGCC plants in general have inherent implementation problems 
due to their integrated nature (Ansolabehere et al., 2007; Maurstad, 
2005). Consequently, several authors have discussed the possi¬ 
bility of integrating the following: the air separation unit (ASU) 
and the gas turbine (GT) streams; and the production of steam 
from several hot stream sources by means of the heat recovery 
steam generation (HRSG) system. Two aspects of ASU-GT inte¬ 
gration are considered: the use of the GT compressor instead of 
an independent ASU compressor to supply the high pressure air 
required in the ASU; and the reinjection of waste N 2 from the 
ASU into the GT combustor for flame temperature control. All 
authors agree that use in the GT combustor of waste N 2 from 
the ASU leads to higher electricity production and less nitrogen 
oxides emissions, and that the use of the GT compressor instead 
of an independent ASU one leads to power consumption sav¬ 
ings (Emun, Gadalla, Majozi, & Boer, 2010; Frey & Zhu, 2006; 
Wang, Qiu, Wu, & Li, 2010). Regarding the HRSG system in par¬ 
ticular, or the heat and power integration of the entire plant in 
general terms, several authors have applied mathematical opti¬ 
misation tools and heat integration heuristics (pinch analysis) to 
provide insights into the required heat exchangers network (HEN) 
(Bassily, 2005; Elia, Baliban, & Floudas, 2010; Gassner & Marechal, 
2009a), as well as the number and kind of pressure levels in the 


steam turbine (ST) cycles. There is no one optimal response for all 
cases. 

In summary, this paper examines the context and the current 
IGCC power plant challenges, namely co-gasification, C0 2 capture 
technologies and integration options. The subject is tackled through 
conceptual modelling. First, we explain the modelling approach and 
provide details of each modelled unit, which are then compared 
with industrial plant data. Second, in Section 3, we describe the 
behaviour of the most relevant plant parameters. The aim here is 
to increase our understanding of the parameters influence on the 
performance of the entire system. We also compare different co¬ 
gasification and co-production scenarios. Finally, Pareto curves are 
introduced as a comprehensive approach that helps the decision¬ 
making task in comparisons of different scenarios. 

2. IGCC power plant modelling 

In this case, the objective of process design is to generate the 
most efficient, economic and environmentally friendly scenario 
for a standard IGCC power plant configuration with and without 
carbon capture units. The model allows for the addition of novel 
units and for changes in unit parameters, as required in the pro¬ 
cess design. In this section, we describe the plant units and their 
modelling assumptions. Then, for calibration and validation pur¬ 
poses, we compare the results of the model with real plant and the 
literature data. The model has been implemented in Aspen Plus® 
commercial software. 

2.1. IGCC with carbon capture units: base case description 

The main steps in an IGCC process with carbon capture units 
can be grouped into five blocks that pursue a common objective: to 
transform a solid fuel into a valuable gas that can be burnt, cooled 
down and expanded in a combined cycle (CC), or transformed and 
sold as pure H 2 . The model was calibrated using data from the Elco- 
gas power plant and from studies by Desideri and Paolucci (1999); 
Hamelinck and Faaij (2002), Chiesa, Consonni, Kreutz, and Williams 
(2005) and Descamps, Bouallou, and Kanniche (2008). Our previous 
work (Perez-Fortes, Bojarski, Velo, & Puigjaner, 2010), which was 
based on kinetic models, was used to adjust other variables. These 
calibration variables can be changed and optimised. The abovemen- 
tioned five blocks are described here and represented in Fig. 1. 

• Feedstock dust preparation. This involves the processes required to 
adapt solid raw material to the gasifier input conditions, including 
drying and grinding the inlet mixture. We do not consider pre¬ 
treatment options like torrefaction or pyrolysis that enhance the 
raw material’s properties, such as its LHV (through densification) 
or its grindability (Gassner & Marechal, 2009b). Therefore, we 
assume that any type of biomass that goes into the system has 
the appropriate characteristics to be pre-processed like coal. 

• Gasification and waste heat boiler (WHB). This block incorporates 
the gasifier itself and the syngas cooling before the cleaning units. 
The model considers a pressurised entrained flow (PRENFLO) 
gasifier, as entrained bed gasifiers are the most frequently used 
reactors in IGCC plants (Highman & van der Burght, 2003). PREN¬ 
FLO gasifiers allow ashes to melt into the slag. Consequently, 
syngas is obtained with notably low emission of fly ashes. More¬ 
over, they achieve a higher char conversion than other types 
of gasifier, which implies that the process is more efficient. 
Entrained flow gasifiers require the fuel to be converted into small 
particles (in the order of microns), which means that a previ¬ 
ous milling step is essential. Moreover, they usually use “pure” 
oxygen (enriched air of around 85% on a molar basis) and steam 
as gasification and moderator agents. The WHB uses waste heat 
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Consumables 



Fig. 1 . IGCC power plant layout, indicating main mass flows and IGCC operations interconnectivity. Inlet and outlet streams are marked in red. Main stream, with final 
products, are marked in bold. The different integration systems are marked with different colors. (For interpretation of the references to color in this figure legend, the reader 
is referred to the web version of the article.) 


from the process of cooling down the syngas before it enters the 
cleaning units to produce steam for the CC. 

• Syngas cleaning. Syngas needs to be cleaned and conditioned 
before it can meet the requirements of its final use, which may 
range from chemical or fuel production to power and/or heat. In 
the IGCC configuration, the objective is to clean the gas before 
it is combusted in the GT combustor to avoid, as far as possible, 
nitrogen and sulphur oxide emissions to the atmosphere as well 
as turbine damages. In the case of H 2 generation, the objective 
is to separate C0 2 from H 2 . The high H 2 content stream is fur¬ 
ther purified and sold as pure hydrogen on the market. In the 
first step of both configurations, solids are removed from the 
gas in a ceramic filter. Secondly, the syngas is placed in a ven¬ 
turi scrubber in contact with water that removes cyanide, halide, 
acid (mainly H 2 S) and basic (mainly NH 3 ) species. Thirdly, before 
its combustion or H 2 purification, the gas goes through an acid 
species removal step that involves a COS hydrolysis reactor and 
an MDEA absorber. The polluted water from the venturi scrubber 
is pre-treated in a sour water stripper before it is transported to a 
waste water treatment plant. The sulphur that is present as H 2 S in 
the syngas is recovered as liquid sulphur in a Claus plant. The syn¬ 
gas cleaning layout follows the configuration of the Elcogas plant 
(Coca, 2003). In the case of the hydrogen generation application, 
the syngas cleaning block includes the carbon capture block. In its 
pre-combustion configuration, it requires a water gas shift (WGS) 
reactor that produces H 2 and C0 2 , followed by a process for C0 2 
capture to separate C0 2 from H 2 . In this simulation, C0 2 removal 
is performed using physical absorption through a Rectisol process 
that employs pure methanol as a solvent. The H 2 is then further 
purified by means of a pressure swing adsorption (PSA) system. 
Finally, the C0 2 stream is sent to a compression system to be 
liquefied for its transport. 

• CC and heat recovery steam generator (HRSG). The CC employs 
Brayton and Rankine cycles to produce electricity. The steam 
cycle uses the waste heat from the system to produce steam at 
different working pressures (high pressure, intermediate pres¬ 
sure and low pressure [HP, IP and LP]). Specifically, the HRSG is 
a HEN that uses waste heat from the flue gas in the gas cycle to 
produce steam. The GT and ST are connected through this flue 
gas, which needs to be cooled down before it is released into the 
atmosphere. Thus, the processes can be integrated. 

• Air separation unit(ASU). The ASU produces oxygen at the required 
purity to be sent to the gasifier as a gasification agent. Pure 
N 2 and waste N 2 are also obtained and are mainly used for 


raw material transportation, inertization, gasification, tempera¬ 
ture moderation, pressurization and nitrogen emission reduction, 
respectively. 

To sum up, the five blocks in the gasification plant are connected 
as follows: the raw material that enters the system is firstly dried 
and crushed. The obtained feedstock dust is then transformed into 
syngas in a PRENFLO gasifier. Before it enters the gas purification 
block, the syngas is cooled down in a WHB that allows the heat to be 
used in the ST. Next, the syngas is cleaned by removing solids. Basic 
and acid trace components are eliminated by means of absorp¬ 
tion cleaning processes. The clean gas is finally sent to a CC. Liquid 
sulphur is obtained as a by-product. 

This type of plant has higher overall efficiency when its underly¬ 
ing blocks are integrated. The integrated factors considered in our 
approach are as follows (see Fig. 1): 

• All the oxygen that the gasifier requires comes from the ASU, 
which uses compressed air from the CC compressor, rather than 
independent compression. 

• After the ASU process, two N 2 streams are obtained: a relatively 
pure one, and waste nitrogen. The pure N 2 stream is used for 
transportation and pressurisation in the feedstock dust prepa¬ 
ration block, and for temperature moderation in the gasifier. In 
turn, the residual N 2 from the ASU is used to control NO* for¬ 
mation in the GT combustion chamber. In plants such as Elcogas, 
the N 2 network is not only used as a service, but also as a work¬ 
ing tool for pneumatic pumps, fly ash discharge, filter cleaning or 
slag transportation. These applications have not been considered 
since they involve relatively small fluxes, which are of no concern 
in a preliminary design approach. 

• In the steam network, the ST cycle provides steam for clean gas 
saturation before combustion in the GT, and for the gasification 
process. It also provides water for the venturi scrubber. In turn, 
heat from syngas cooling, in the WHB and the HRSG is used to 
produce steam in the CC. 

Fig. 1 shows the overall mass balance of the plant. The main 
consumers of the net power production are GT cycle compression, 
the ASU operation, C0 2 compression in the Rectisol process and C0 2 
compression for liquefaction. The main technical data considered in 
the model are summarised in Table 1. In the following subsections 
we use the modelling strategy to provide more details for each one 
of the considered units. 
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Table 1 

Technical data for the IGCC plant model. 


Operating conditions 


Gasification 

Feedstock 

2600 ton/day 

Crusher fineness 

50-60 |jim 

Feedstock MC 

2 % 

P of lock hoppers 

30 bar 

Tgasif 

1400-1500°C 

Pgasif 

25 bar 

0 2 / 0 2 stoich (mole basis) 

Raw material dependant 

H 2 O/Q n (mass basis) 

Raw material dependant 

T before the WHB 

850°C 

T before gas cleaning 

235 °C 

Max. P before CC 

3 bar 

Combined cycle (CC) 

Air ^ 

Adjusted to regulate T gasif and TOT 

HP/IP/LP 

127/35/6.5 bar 

HP/IP/LP T before ST 

510/520/265 °C 

Compressor pressure ratio 

15.7 

T cg before combustor 

300 °C 

TOT 

540 °C 

P of the saturation column 

20.8 bar 

T of the flue gas 

100°C 

T of the air to ASU 

130°C 

Consumables 

NaOH (15 mass%) 

Syngas and cleaning water dependant 

H 2 SO 4 (96mass%) 

Cleaning water dependant 

Natural gas 

Feedstock dependant 

h 2 o 

Syngas and raw material dependant 

Air separation unit (ASU) 

Oxygen purity 

85 vol% basis 

Radiation losses 

2 % adiabatic reactor 

To 

15 °C 

Po 

0.93 bar 


2.1.1. Simplifying hypotheses and general assumptions 

As mentioned above, we do not take all the steam and N 2 net¬ 
works in real power plants into consideration in this study. Instead, 
we focus on the main streams. The IGCC modelling approach does 
not consider the slag cooling and collection system or an exhaus¬ 
tive model of the cooling tower. The design of the waste water 
treatment plant is also outside of the scope of this paper. The Elco- 
gas power plant employs gasifier technology from Krupp-Koppers 
and a CC from Siemens. Their output data, as well as input and 
output data from their gas cleaning units with no carbon capture, 
were used to calibrate the model. As a general physical property 
method to calculate the thermodynamic and transport properties 
of streams, we opted for the Peng-Robinson equation of state with 
the Boston-Mathias alpha function (PR-BM) method. This equation 
of state is recommended to model gas phase systems at medium 
to high pressure. NBS/NRC Steam Tables were used for the ST 
with water as a working fluid, since they are implemented in the 
same way as any other equation of state (AspenTech, 2010). The 
ELECNRTL method was chosen to deal with the electrolyte sys¬ 
tems, namely the venturi scrubber, the sour water stripper and 
the MDEA absorber. The NRTL method was used for the physical 
absorption process, which is the Rectisol process in this study, fol¬ 
lowing the VL equilibrium through Henry’s law as in the ELECNRTL 
property package, but with the difference that no species ionisation 
occurs. These last models have been replaced by meta-models (in 
this case, component splitter models) which are calculated from 
the exhaustive models, due to the excessive time required for 
their resolution in sensitivity analyses or variables optimisation 
calculations. The LHVs were estimated using Aspen Plus® prop¬ 
erty sets (QVALNET). The stoichiometric value of air was calculated 
using the COMB-02 property set from Aspen Plus® (AspenTech, 
2010 ). 


Table 2 

Ultimate and proximate analyses of raw materials (Coca, 2003; ECN-Biomass, 2010). 


% Mass basis (dry) 

Coal 

Petcoke 

Orujillo 

Ultimate analysis 

C 

41.07 

88.40 

50.00 

H 

2.81 

3.34 

6.50 

O 

7.51 

0.02 

36.30 

N 

0.92 

2.04 

0.80 

S 

1.05 

5.91 

0.10 

Cl 

0.04 

0.00 

0.20 

Proximate analysis 

Moisture 

11.80 

7.00 

7.60 

Ashes 

46.60 

0.28 

6.30 

Fixed carbon 

32.05 

85.74 

21.30 

Volatiles 

21.35 

13.98 

72.4 

HHVdry (MJ/kg) 

13.58 

32.65 

20.38 


2.1.2. Feedstock dust preparation 

This block includes feedstock dust preparation, drying, pres- 
surisation and gasifier feeding. The solid inlet stream is defined 
in Aspen Plus® as the addition of non-conventional components. 
The composition of the base case coincides with the base com¬ 
position of Elcogas for validation purposes. Indeed, it corresponds 
to a 50:50 mixture of coal and petcoke on a mass basis. The coal 
comes from the local Encasur mines and is sub-bituminous with 
high ash content. The petroleum coke has high sulphur content and 
comes from the Puertollano Repsol refinery, where it is obtained as 
a by-product. The residual biomass used as feasible waste for the 
co-gasification analyses in the IGCC power plant is olive pomace or 
orujillo, which is the solid residue that remains after pressing olives. 
Table 2 shows the main data for the feedstock composition. The 
higher heating value (HHV) of the three fuels is taken into account 
when the data are introduced into the simulator. In addition, lime¬ 
stone is fed into the gasifier as an additive to reduce the ash fusion 
temperature. The limestone contains 95% CaC0 3 and 5% ashes. It is 
roughly 2-3% in weight of the total feedstock introduced into the 
system, which we consider to be fixed regardless of the blend used. 

The model allows for mass composition changes in the coal, pet¬ 
coke, and orujillo blend using a Fortran code (introduced by means 
of a calculator block) to calculate the ultimate and proximate anal¬ 
yses and the HHV of the blend, on the basis of each fuel flow. The 
maximum feedstock flow is 2600tons/day, which corresponds to 
100% of the gasifier/plant load. The Fortran code also controls the 
plant load. The dust preparation itself takes place in a crusher unit 
so that the inlet fineness can be changed from about 100 mm to 
50-60 |jim. Additional information to be introduced is the mixture’s 
grindability, in terms of the Bond work index. This is calculated in 
a calculator block using a Fortran code and taking into account the 
mass proportion of each component. The Bond index values are 
73.80 kWh/ton for petcoke and 11.37 kWh/ton for coal. The oru¬ 
jillo is taken to have the same value as the coal. Since limestone is 
not in the consulted database, we selected the Bond work index for 
dolomite, which is 11.31 kWh/ton (Perry & Green, 1997). The fuel 
drying step involves two reactors: a combustion chamber, in which 
hot gases are generated, and the drier itself, in which the feed¬ 
stock stream is dried by solid-gas contact. A Gibbs reactor model 
is used for the combustion reactor, where natural gas is burnt with 
an excess of air. The excess of air is established using a calculator 
block. The air needed for the feedstock drying has been set at 1.2 
times the required stoichiometric oxygen. This reactor operates at 
0.93 bar and the temperature is estimated as a fraction of the flame 
temperature. The drier is modelled using a stoichiometric reactor 
with the feed stream as the limiting reactant, to achieve a final 
moisture content of 2% on a mass basis, as stipulated in a calculator 
block. Before the powdered fuel enters the gasifier, it is separated 
from the inert gas in bag filters. These are simulated with a flash 
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separator. Before the feed enters the gasifier, it is pressurised at 
30 bar inside a lock hopper system. This system is simulated in a 
simplified way as a mixer, taking into consideration the pure N 2 
stream. Finally, the mixture is conducted to the gasifier. 

2.1.3. Gasification and the WHB 

Gasifier modelling approaches were reviewed in Villanueva, 
Gomez-Barea, Revuelta, Campoy and Ollero (2008) and in Perez- 
Fortes et al. (2009). They range from the macroscopic level that 
includes thermochemical equilibrium approaches, to the micro¬ 
scopic level, which incorporates chemical kinetics. The most 
extended model for an entrained bed gasifier is a Gibbs reactor, 
as discussed by Emun et al. (2010) or Kunze and Spliethoff (2010). 
Therefore, we use a Gibbs reactor in this study. Gasification reac¬ 
tions can be described by many different equations. The most 
characteristic gasification reactions are given in Eqs. (1)-(10), and 
correspond to reactions with oxygen, steam, C0 2 and hydrogen. 
Although the raw material composition is based on C, FI, O, N, S and 
Cl, the main reactions only involve C as reactant. Eqs. (l)-(3) cor¬ 
respond to combustion, Eq. (4) is the Boudouard reaction, Eq. (5) is 
the primary WGS and Eq. (6) is the WGS. Eq. (7) is the WGS reaction 
with gas species only, and Eq. (8) is the methanation reaction. Eqs. 
(9) and (10) correspond to steam and dry reforming. As reactions 
with oxygen are complete under gasification conditions, the main 
role can be seen in Eqs. (4)-(8) in terms of equilibrium (Fermoso 
et al., 2010; Highman & van der Burght, 2003). 


C + 0.5O 2 CO (1) 

CO + 0.5O 2 C0 2 (2) 

H 2 + 0.5O 2 ^H 2 O (3) 

C + C0 2 2CO (4) 

C + H 2 0*>C0 + H 2 (5) 

C + 2H 2 0^C0 2 +2H 2 (6) 

CO + H 2 0 C0 2 + H 2 (7) 

C + 2H 2 ^CH 4 (8) 

CH 4 + H 2 0 CO + 3H 2 (9) 

CH 4 + C0 2 2CO + 2H 2 (10) 


The main gasifier operation conditions, the pressure and tem¬ 
perature (Pgasif and T gasi f ), are driven by the GT working pressure 
and ash melting point, respectively. Thus, the GT pressure sets the 
gasifier pressure at 25 bar. The gasification temperature is speci¬ 
fied by the requirement for ashes to be separated as molten slag. 
According to Song et al. (2009), this temperature should be between 
1400° C and 1500° C for coal, depending on the limestone con¬ 
tent. In the model, we have established this temperature interval 
as a condition in gasification when the raw material composition is 
changed. Specifically, a temperature of 1450° C is fixed as a design 
specification in Aspen Plus®. 

The feedstock from the feedstock dust preparation step is a non- 
conventional stream. Consequently, before it is introduced into the 
Gibbs reactor it should be divided into its elements, since there 
is no Gibbs energy information for a non-conventional component. 
Therefore, before the gasifier itself, a yield reactor is used to decom¬ 
pose the raw material into C graphite , H 2 , 0 2 , N 2 , S, Cl 2 , water and 
ashes, according to its ultimate analysis and the moisture content 
from the proximate analysis. T gasi f is estimated by introducing into 
the Gibbs reactor the heat of reaction associated with the feedstock 
stream decomposition and the heat released into a component 
separator that separates the raw syngas from the slag (ashes and 
limestone). 


Oxygen from the ASU, which is the main gasifier agent, and IP 
steam and N 2 , which act as moderators, are also introduced into the 
Gibbs reactor. Oxygen and steam flowrates are estimated by spe¬ 
cific ratios that depend on the raw material composition, according 
to its ultimate analysis characteristics. The oxygen ratio is defined 
as the ratio between the oxygen flowrate and the raw material 
mixture’s stoichiometric oxygen, after decomposition into its com¬ 
ponents on a molar basis. This quotient is known as the equivalence 
ratio (ER). The steam ratio is defined on a mass basis and has the 
mass flow of Cg raphite as a denominator. The values of these ratios 
are 0.42 and 0.16, which correspond to the base case (50% of coal 
and 50% of petcoke on a mass basis). They have been estimated by 
adjusting the gasification temperature to 1450° C and using the 
Elcogas input conditions, respectively. 

The WFIB is introduced after the production of syngas. This 
HEN uses the heat from cooling the syngas before it enters the gas 
cleaning units. Syngas cooling takes place in two steps: firstly, the 
temperature of gasification gases is reduced when they are mixed 
with the quench gas, which is a fraction of the cooled syngas itself 
that is recycled for this purpose (at around 235 °C). In this step, the 
gases are cooled to 850° C which is close to the heat exchanger’s 
material temperature limits (Coca, 2003). Secondly, in the WHB, 
the gas heat is recovered in a HP heat exchanger, so that the gas 
is cooled to 380° C and HP steam is generated. Finally, the gas is 
moved to an IP heat exchanger that cools it to 235° C and gen¬ 
erates IP steam. After this last step, the quench gas is recycled to 
be mixed at the gasifier outlet. The recycled mass flow is deter¬ 
mined by a design specification block, in which the temperature of 
850 0 C is attained by modifying the quench gas flowrate. Feedstock 
conditioning, gasification and syngas cooling are schematically rep¬ 
resented in Fig. 2. 

2.1.4. Syngas cleaning 

The aim of this block is to obtain a clean gas before combustion 
in a CC to avoid pollutant emissions, as far as possible. In addi¬ 
tion, C0 2 may be separated from H 2 to further purify the hydrogen 
stream and thus produce pure hydrogen to be sold to the market. 
The basic cleaning process for both final applications is made up of 
cleaning units that work at HP (over 20 bar). The syngas is cleaned 
in the venturi scrubber and the COS hydrolyser-MDEA absorber. 
By-product streams are cleaned in the sour water stripper and in 
the Claus plant. The layout of the cleaning units is represented in 
Fig. 3. 

• Venturi scrubber. This unit is placed after the gasifier and after the 
removal of solids. Its main objective is to abate ammonia, cyanide 
and halide species. Some H 2 S is also absorbed, although this is not 
the main purpose of the scrubber. This process unit allows contact 
between water and syngas. Thus, species transfer occurs accord¬ 
ing to the VL equilibrium, following Henry’s law and considering 
the electrolyte reactions that take place in the liquid phase. The 
parameters are reported elsewhere (Perez-Fortes et al., 2010) and 
can be found in the Aspen Plus® Properties data bank. The ven¬ 
turi scrubber is modelled using a RadFrac column with a constant 
pressure of 23.6 bar. It uses a 15% solution of NaOH (by weight) 
that acts as an acid capturer. Water and NaOH consumption are 
proportional to the syngas flowrate that enters the venturi scrub¬ 
ber: 9% and 0.12% of the raw gas flowrate have been determined 
respectively. The proportionality factors are introduced into the 
model through two Fortran codes, by means of two calculator 
blocks. Polluted water goes to a sour water stripping system to 
be pre-treated before the final water disposal process. This sys¬ 
tem works at nearly atmospheric pressure (1.5 bar). Hence, water 
is depressurised in a flash vessel before it enters the first column. 
The flash vessel has been modelled as a two-phase flash reactor 
that works at 1.5 bar and 53 °C. The clean syngas that is obtained 
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Fig. 2. Feedstock conditioning (dust preparation and drying) and gasification steps. 


after the venturi scrubber is cooled to around 140° C before the 
COS hydrolyser. 

• Sour water stripping system. The main purpose of this system is 
to pre-treat the used water in the venturi scrubber by desorb¬ 
ing NH 3 , H 2 S, C0 2 , HCN and HC1. It consists mainly of two units: 
an acid stripper that abates acid species and a basic stripper that 
abates basic species. Again, this system is modelled using elec¬ 
trolytes and considering the VL equilibrium by applying Henry’s 
law. The desorbing columns are pH regulated by adding acid and 


basic solutions: H 2 S0 4 and NaOH with a purity of 96% and 15% by 
weight, respectively. The solution flowrates depend on the water 
flowrate from the venturi scrubber. The relationships are intro¬ 
duced into the model through Fortran calculator blocks: the acid 
solution is 0.3% of the overall water mass flowrate that needs to 
be treated. Similarly, the basic solution corresponds to 1.2% of the 
mass flowrate that needs to be treated. The strippers are mod¬ 
elled with RadFrac columns with condenser and reboiler units. 
After the stripping process, a sour gas stream is obtained that is 



CLEAN GAS 
TO CC 


CLEAN GAS 
TO PRODUCE 
H 2 


Fig. 3. Cleaning units layout indicating main mass flows. Electricity and H 2 production superstructure is also represented. Note that the inlet flowrates are marked in red. 
Superstructure choices used in this work are marked with dotted and green lines. (For interpretation of the references to color in this figure legend, the reader is referred to 
the web version of the article.) 
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directed to the Claus plant for sulphur recovery. The ionic dis¬ 
sociation reactions that are considered in both columns and in 
the venturi scrubber are described elsewhere (Perez-Fortes et al M 
2010 ). 

COS hydrolyser and MDEA absorber. The purpose of these systems 
is to remove syngas sulphur compounds. The COS transformation 
into H 2 S has been modelled using a stoichiometric reactor with 
COS fractional conversion of 0.99. This ratio is based on a kinetic 
model described in Perez-Fortes et al. (2009). The MDEA-water 
solution acts as a chemical absorber for FI 2 S at TIP and low tem¬ 
perature (22 bar and 33 °C). Absorption and desorption columns 
are simulated, with the subsequent MDEA lean stream recircula¬ 
tion. The MDEA solvent is a 50% water solution on a mass basis. 
The amount of MDEA solution that is recycled is calculated by 
a sensitivity analysis for the base case, using Elcogas data. Due 
to the MDEA losses, the net consumption of MDEA is around 
3 kg/h. The absorption-desorption columns are modelled using 
RadFrac models. The desorber works at 1.5 bar and relies on a con¬ 
denser and a reboiler. The MDEA system considers the presence of 
electrolytes. The ionic dissociation reactions and the equilibrium 
expressions are described elsewhere (Perez-Fortes et al., 2010). 
Claus plant. The aim of this plant is to recover sulphur through 
H 2 S conversion into elemental sulphur. A secondary aim is to 
convert NEI 3 and EICN into N 2 . The Claus plant works with two 
polluted streams: Claus gas from the MDEA absorber and sour 
gas from the sour water stripping process. Both gas streams are 
at 1.5 bar, and the process takes place at 0.93 bar. The proce¬ 
dure relies on two parallel kilns and two series reactors that 
are catalysed by alumina. The Claus plant process finishes with 
a catalysed hydrogenator that is used to increase the overall sul¬ 
phur removal efficiency through residual gas recycling. Thus, 
it provides another opportunity to remove the remaining sul¬ 
phur from the gas. The first thermal and catalytic stages have 
been modelled with stoichiometric reactors, based on the results 
obtained in a previous kinetic-based model (Perez-Fortes et al., 
2010). The thermal stage takes place in one single reactor, which 
is represented by Eqs. (11)-(14). The catalytic stage involves Eq. 
(14), at two different temperatures: 270° C and 210 °C. There are 
two streams for air in the Claus kiln; 33% of the stoichiometric 
oxygen that is needed to combust the H 2 S present in the gas; and 
the stoichiometric air that is required to combust the rest of the 
gas components. 


H 2 S + 1.50 2 -3 

► so 2 + h 2 o 

(ii) 

NH 3 -> 0.5N 2 

T 1.5H 2 

(12) 

HCN + H 2 0 -> 

1.5H 2 + 0.5N 2 + CO 

(13) 

2H 2 S + S0 2 -* 

3S + 2H 2 0 

(14) 


The air that is fed to the Claus kiln is calculated as the fraction 
of the stoichiometric air that is required to burn the quantity of 
H 2 S present in the gas mixture, which is set in the simulation 
by means of a calculator block. The assumed degrees of advance 
of the different chemical reactions are: total conversion of the 
0 2 from Eq. (11), and 0.96 of NH 3 and HCN conversions in Eqs. 
(12) and (13). In the case of Eq. (14) the EI 2 conversion is con¬ 
sidered to be 0.55, while in the case of the two catalytic stages, 
a conversion of 0.8 and 0.6 is proposed for the S0 2 . Finally, the 
hydrogenator consumes a small portion of clean syngas (0.5%) to 
transform unrecovered S and untransformed S0 2 into H 2 S, which 
then returns to the gas cleaning process. The hydrogenator has 
been simulated with a Gibbs reactor. After each thermal and cat¬ 
alytic stage, S is recovered by means of flash vessels that cool the 
stream flow to collect the liquid sulphur that is produced. 

WGS and C0 2 capture. The pre-combustion configuration was cho¬ 
sen for removing C0 2 to obtain H 2 . 


- The first step before C0 2 capture is the conversion of CO into 
H 2 and C0 2 . This is carried out in the WGS reactors. This step 
is modelled following the configuration described in Chiesa 
et al. (2005). The shifting transformation takes place in two 
adiabatic reactors with Eq. 7 in equilibrium at two different 
temperatures: 400° C and 210°C. Downstream from the first 
reactor, the gas stream is cooled to meet the conditions of the 
second reactor. Downstream from the second reactor, the gas 
stream is cooled from 320° C to 200° C so that it can meet 
the C0 2 removal unit conditions. Note that the WGS reaction 
is exothermic. The first reactor consumes IP steam to reach the 
high temperatures that are required. The amount of steam is 
fixed through the j3 parameter, which is defined as the steam- 
to-CO ratio and is set to a value of 1.25, according to Descamps 
et al. (2008). This step is modelled in Aspen Plus® with two 
equilibrium reactors. 

- C0 2 removal takes place using a Rectisol process, with 
methanol as the absorption solvent. The solvent expansion 
between columns takes place by means of three pressure drums 
(modelled as three flash vessels in Aspen Plus®), which recover 
the C0 2 and allow for solvent recycling. These three units have 
been established to obtain the desired methanol purity. The 
high H 2 concentration stream is obtained at the top of the 
absorption column, where the purity is 80% on a molar basis. 
The final C0 2 stream comes from the decompression step and 
from the top of the desorption column. The purity is 70% on 
a molar basis. The final pressure of the EI 2 stream is 21.7 bar, 
while that of C0 2 is lObar. The absorption column works at 
low temperatures of around -30°C. The desorption column 
works at 1.2 bar (note that the final C0 2 stream is compressed 
to the abovementioned lObar), and includes a condenser and 
a reboiler. The amount of methanol required is proportional to 
the mass flowrate of C0 2 contained in the WGS product stream: 
it is 4.5 times the C0 2 mass flowrate that needs to be treated. 
This value is introduced into the model through a calculator 
block. The variable values of the Rectisol unit are based on the 
work of Elamelinck and Faaij (2002) and Descamps et al. (2008). 
The two main final streams contain other species, the H 2 stream 
has a fraction of N 2 (around 15% on a molar basis), and the C0 2 
stream has a small fraction of syngas compounds. The methanol 
flowrate and the number of plates in the RadFrac models are 
fixed, but the N 2 stream can limit the Rectisol process. This is 
because the partial pressure of the C0 2 stream can be limited 
by an increase in the N 2 fraction, which is not large enough to 
displace the VL equilibrium and be absorbed. On the H 2 side, 
the presence of more N 2 , which flows together with the EI 2 
stream, reduces the stream purity and thus means that further 
treatment is required in the PSA system. 

- After the Rectisol unit, the EI 2 can either be further purified 
or used to feed an ideal CC that works with hydrogen as a com¬ 
bustible. The stream should be heated in both applications: up 
to 40° C to meet the PSA requirements, or up to 150° C to 
meet the syngas conditions that are required after the MDEA 
absorber and before the entry to the CC. The configuration and 
parameter values for the PSA system for further purification of 
the H 2 stream are from a study by Elamelinck and Faaij (2002). 
The system is based on two reactors that use two different solid 
beds: activated carbon and a zeolite molecular sieve. The trade¬ 
off in these systems is the number of reactors vs. the flue gas 
recycled stream to enhance H 2 recovery. After this process, the 
flue gas can be released to the atmosphere or compressed and 
used in the CC (see Fig. 3, PSA adsorption block). Its LEIV is 
around 8 MJ/kg. Each PSA unit is modelled with a component 
separator: in the first one it is assumed that all the C0 2 and all 
the H 2 0 are released from the main H 2 stream, and in the sec¬ 
ond one the remaining EI 2 corresponds to 84% of the total H 2 
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Fig. 4. Simplified flowsheet of the GT cycle. 



Fig. 5. Simplified flowsheet of a ST cycle. 


flow. The final H 2 pressure is around 0.95 bar, with a purity of 
99.99% on a molar basis. 

- The C0 2 stream needs a liquefaction step to ease trans¬ 
portation and for geological storage. The liquefaction system 
follows Desideri and Paolucci (1999). The operating pressure is 
assumed to be 150 bar. The compressor has an isentropic effi¬ 
ciency of 82% and an outlet temperature of35 0 C before pipeline 
transportation. Before its compression, the water fraction of the 
C0 2 stream is separated from the main stream with a flash ves¬ 
sel to avoid problems in the pipeline. This is why the lObar 
stream is cooled from 183 0 C to 25 0 C after the desorption col¬ 
umn from the Rectisol system. The final C0 2 stream has a purity 
of 95.5% on a molar basis. 

2.1.5. CC and HRSG 

Regarding the overall CC modelling strategy, we follow the 
approach found in Zhu (2004) and Ongiro, Ugursal, Altaweel, 
and Blamire (1995). The GT cycle is composed of a compres¬ 
sor/combustion chamber/turbine system that uses clean gas. A 
saturation column situated before the clean gas combustion pro¬ 
cess saturates the stream with vapour and adds nitrogen. A 
simplified flowsheet of the process is shown in Fig. 4. Turbine and 
compressor stages are modelled with turbine and compressor units 
from the Aspen Plus® library. Several mechanical and operating 
limitations have been disregarded such as choking and surging in 
the compressor (Giampaolo, 2006) or the gas turbine inlet tem¬ 
perature (TIT), which has been considered free and not subjected 
to design limitations. The saturation column is a two-phase flash 
model in which water is added until the relative syngas humid¬ 
ity reaches 100% (Coca, 2003). This is stipulated using a design 
specification. Turbine air cooling has been modelled by taking into 
consideration a 4-stage GT. Thus, the compressor has also been 
modelled as a 4-step process. Each corresponding stage (i) has the 
same pressure loss or gain ratio according to Eq. (15). This is intro¬ 
duced into the model by Fortran code using a calculator block, in 
which PR is the pressure ratio of the turbine or compressor, P 0 is 
the inlet pressure, and n s t is the number of stages. 

P i(<ss = P 0 • PR('/ n «)% (15) 

The combustor is modelled by a Gibbs reactor. Turbine outlet 
temperature (TOT) is a critical value since it influences the amount 
of heat recovered in the HRSG. It is controlled by the air mole 
flowrate of entry into the combustor by using a design specifica¬ 
tion. As in the gasifier, the design specification sets the air mole rate 
as the proportional stream of the stoichiometric air that is needed 
to burn the clean gas mixture and obtain the desired TOT, which 
is set at 540 °C. In our base case, the proportional factor value is 
2.45. However, this is modified by the syngas composition and the 
flowrate variations. The temperature in the combustor is 98% of 
the adiabatic one. Cooling air fractions are established by optimis¬ 
ing the GT block in the range of 10-20% the inlet gas per stage. The 


HRSG is a very complex HEN that considers various stream temper¬ 
atures, and combines and intercalates heat exchangers to produce 
steam at different temperatures and degrees of superheating. A 
simpler approach has been considered in this conceptual design. 
The heat transfer area and heat transfer coefficients are disregarded 
when we calculate the heat produced from the flue gas temperature 
changes. Nevertheless, we consider that the outlet gas tempera¬ 
ture must be higher than the steam temperature that we wish to 
produce, to prevent temperature crossing. The HRSG takes advan¬ 
tage of the syngas temperature drop from 540° C to 100°C. Thus, 
four heat exchangers have been considered to generate HP, IP and 
LP steam and to condensate water which then returns to the feed 
water tank. To maximise the power generated in the ST, the out¬ 
let temperatures for the HP, IP and LP steam sides heat exchangers 
have been estimated and set at 300 °C, 180 0 C and 170 °C. 

Three cycles are considered in the ST cycle. They are composed 
of a water pump, economiser, evaporator, superheater and turbine, 
using analogous units from Aspen Plus® (see a general flowsheet 
in Fig. 5). Water is available at different pressures and it flows from 
its corresponding cycle (HP or IP) down to the LP cycle. There¬ 
fore, the LP turbine produces the highest amount of power, due 
to the large amount of steam that is depressurised. Water mass 
flows for each pressure cycle are calculated by considering the total 
heat that is recovered from the WHB and the HRSG, as well as the 
final estimated temperature before each turbine decompression. A 
calculator block is used for these values assignments, see Table 1 
for more details. Water consumption has to be considered, as net 
consumption varies in the plant and includes the steam fed to the 
gasifier, the steam consumed in the venturi scrubber to clean the 
syngas, the steam flow used in the saturator to saturate the clean 
gas before its combustion in the GT, and the water consumption 
needed in the WGS unit to produce H 2 and C0 2 from CO. The con¬ 
sumption of these streams penalises the heat flow from WHB and 
HRSG to the ST cycle. Finally, the net CC power is obtained by con¬ 
sidering the gross power resulting from the ST and the GT (taking 
into account the air compressor consumption) and subtracting the 
power consumption of the ASU, the crusher, and the compressors 
and pumps in the syngas cleaning unit blocks, including the power 
of the C0 2 capture system when co-production is considered. 

2.1.6. ASU 

This unit produces enriched air by separating oxygen and nitro¬ 
gen to meet the requirements of the gasification and CC blocks. In 
this case, the unit provides an oxygen stream of 85% purity on a 
molar basis. The ASU product streams are enriched oxygen, pure 
nitrogen and waste nitrogen. These streams are pressurised exter¬ 
nally to adapt them to the plant conditions, since the pressures in 
the ASU process are governed by the cryogenic and distillation con¬ 
ditions of the air. We can distinguish several steps in this procedure, 
which have been modelled in a simplified way with a component 
separator, compressors and heating units: 
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Summary of all modelled units. 


M. Perez-Fortes et al. / Computers and Chemical Engineering 35 (2011) 1501-1520 


1509 


Unit 

Aspen Plus® model 

Feedstock dust 

Dust preparation: crusher 

preparation 

Combustion chamber: Gibbs reactor 

Drier: stoichiometric reactor 

Bag filters: two phase (VL) flash separator 

Lock hopper: mixer 

Gasification and WHB 

Gasifier: yield and Gibbs reactors 

Solid removal: component splitter 

WHB: heating units and a compressor for gas 
recycling 

Syngas cleaning 

Venturi scrubber: RadFrac column, two phase (VL) 
flash separator 

Sour water stripper: RadFrac columns 

COS hydrolyser: stoichiometric reactor 

MDEA absorber: RadFrac columns 

Claus plant: stoichiometric and Gibbs reactors and 
two phase (VL) flash separators 

WGS reactor: two Equilibrium reactors 

Rectisol process: RadFrac columns 

PSA unit: component splitter 

Liquefactor: flash separator and compressor 

CC and HRSG 

Turbines, compressors and heating units 

Combustor: Gibbs reactor 

Saturator: two phase (VL) flash separator 

HRSG: heating units 

ASU 

Compressors and heating units, and component 
splitter 


• Air pre-cooling. In this step, air is taken to the purification unit 
temperature. Generally, this is carried out by a refrigeration fluid. 

• Air purification and distillation. The main objective in the purifi¬ 
cation step is to remove impurities, such as water, from the air 
inlet stream. This step is based on solid-gas absorption. As the 
pressure is determined by the GT compressor, the air pre-cooler 
is crucial to obtain the desired efficiency. Then, the distillation 
stage is a cryogenic process in which nitrogen and oxygen can 
be separated through a distillation column. This occurs at about 
13 bar and -165 °C. Oxygen is released as a liquid and nitrogen 
as a gas. The purity of waste N 2 is approximately 98% on a molar 
basis, and that of pure N 2 is about 99.9% on a molar basis. Instead 
of a complex model, a component separator block is used in our 
approach. 

• Gas supply. The outlet streams need to be adapted to the desired 
conditions of pressure and temperature for the gasifier and the 
CC. These are 30 bar for the 0 2 ,20 bar for the waste N 2 , and 50 bar 
for the IP N 2 . The isentropic efficiencies of the compressors are 
assumed to be 0.72. 

For a summary of the modelling strategy described in Section 

2.1, see all the modelled units and their representation in Aspen 
Plus® in Table 3. 

2.2. Model validation 

The model’s performance is compared with available plant data 
to validate its results. This validation step considers the gasifica¬ 
tion and the gas cleaning blocks separately, as well as the plant as a 
whole. The PRENFLO gasifier produces syngas that mainly consists 
of 21% H 2 , 60% CO, 4% C0 2 ,3.5% H 2 0,1% H 2 S, 10% N 2 and COS, on a 
molar basis in the above-mentioned conditions. The modelled syn¬ 
gas flow contains 20% H 2 ,58% CO, 4.5% C0 2 ,5% H 2 0,1% H 2 S, 11% N 2 
and COS, on a molar basis (Coca, 2003). For the gas cleaning process 
with no H 2 purification, Fig. 6 compares industrial information on 
the clean syngas composition and the composition of the Claus gas 
with the composition obtained from the simulation. Both streams 
are the outlet streams of the MDEA absorption-desorption pro¬ 
cess, which is the last cleaning step applied to the main gas stream 



Fig. 6. Comparison of components of clean gas and Claus gas between values calcu¬ 
lated (predicted) by the model and the actually measured data from industry. The 
error bar corresponds to 10% deviation. 


Table 4 

Overall IGCC plant validation data. 



Measured 

Model 

Error % 

Raw gas flow(t/h) 

202.5 

192.3 

-5.05 

LHV clean gas (MJ/kg) 

9.911 

9.547 

-3.68 

Gross power (MW) 

318 

315 

0.94 

ASU consumption (MW) 

29.0 

27.4 

-5.52 

Net power (MW) 

283 

285 

0.71 

Sulphur (t/h) 

3.11 

3.32 

6.59 


before the CC. The measured and predicted values are on the x and 
y axes. Therefore, the points that are represented reflect the com¬ 
parison between these values. The points closest to the quadrant 
division line are the most precise. If we consider the IGCC power 
plant with no carbon capture units, we can also verify the accuracy 
of the model in terms of the net power, byproducts and clean gas 
characteristics, as shown in Table 4. 

For the feedstock dust preparation, gasification and clean gas 
blocks together, the clean gas and Claus gas composition compar¬ 
ison is within an error margin of 10%. Moreover, the comparison 
results in Table 4 indicate that the biggest discrepancy is in the 
estimation of the sulphur recovered, the flowrate of gas produced 
and the ASU consumption. We can deduce that the underestimation 
of the LFIV of the clean gas has been compensated by a higher effi¬ 
ciency in the CC and lower auxiliary consumptions, rather than the 
ASU one. The PSA unit is directly validated by comparing the model 
results with the data from Chiesa et al. (2005) and Descamps et al. 
(2008). Concerning the WGS reactors, a value of 1.25 is used. The 
CO conversion rate is 91%, in contrast to the 92% rate found in the 
study by Descamps et al. (2008). Concerning the Rectisol system, 
a C0 2 absorption rate of 95% is referred to in the literature, while 
a value of 99% is obtained in the simulation due to the modelling 
assumptions. All in all, the results show that in a conceptual stage 
this modelling approach represents an accurate enough approxi¬ 
mation of an IGCC power plant with C0 2 capture, with a maximum 
error of 10%. This is well within the error range of a preliminary 
design stage, which stands at 15-20% (Wells and Rose, 1986). 

2.3. Metrics 

To deal with the varying performance of the results, we consider 
different efficiency parameters according to the final product of the 
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IGCC plant. The parameters are defined as follows: 


• Cold gas efficiency ( CGE ) (Uson, Valero, Correas, and Martinez 
(2004) is defined as the ratio between the chemical energy con¬ 
tained in the raw gas (rg) (which is what the syngas is called when 
it comes out from the gasifier) and the total chemical energy con¬ 
tained in the fuel. Both of these values are measured using the LHV 
(Eq. 16). In the case of fuel blends, the LHVf ue i is calculated from 
the contribution of each i th fuel on a mass basis (Wj), as in Eq. 


(17). 



CGE = 

TTlrgLHVrg D/ 

(16) 


r u\ / 

mfeed ' LHV feed 

LHV fuel 

= WjLHVj 

(17) 


• The energy efficiency of the combined cycle ( Effcc ) is defined as 
the ratio between the net power of the CC (considering the energy 
consumed by the compressor and the cleaning units) and the total 
chemical energy contained in the clean gas (eg) which is fed to the 
gas turbine, see Eq. (18). According to the co-production options, 
this clean gas can be syngas, hydrogen or PSA purge gas. 


Effcc = 


Power net 
pi eg ■ LHV C g 


(18) 


• By considering the plant as a whole, the total profitable energy 
that enters the system is related with the final electricity that 
goes to the grid. This energy efficiency ( Effglobal ) is calculated as 
in Eq. (19). 


Eff E 


global : 


Power ne t % 
mfeed ' EHV/ ee d 


(19) 


• Similarly, in the case of H 2 production, the global efficiency is 
calculated as in Eq. (21). In the same way as in the CC, the partial 
efficiency of the carbon capture process can be calculated as in 
Eq. (20). 


EffcCapture — 


mu 2 ■ EHV h 2 d/ 
m cg ■ LHV cg ° 


Effglobal^ 2 ~ 


mu 2 • LHV h 2 ^ 
mf ee d • EHVf ee d 


( 20 ) 

( 21 ) 


• Eqs. (16) and (18) evaluate the performance of the gasifier and of 
the CC separately, while Eff g i oba i provides a metric for the overall 
installation. To relate the former definitions, the efficiency asso¬ 
ciated with the gas cleaning operations ( Eff d ) can be defined as 
in Eq. (22), 


Effd = 


m rg ■ LHV rg 

meg • LHVcg 


( 22 ) 


Consequently, Eff g i oba i can be calculated from the partial effi¬ 
ciencies as in Eq. (23). 


Effgiobai = Effcc ■ Eff el • CGE (23) 

The same reasoning can be applied to the H 2 production, result¬ 
ing in Eq. (24) 


Eff global^ = Eff^Capture ■ EJf d ■ CGE (24) 

• In the case that syngas is used for co-production, namely the 
production of hydrogen and electricity at the same time, the effi¬ 
ciency calculation should combine the efficiencies related to both 
electricity ( Eff globa i ) and hydrogen {Eff g i obal ^ ) production. Thus, 
the final global efficiency, called the total efficiency, is defined as 
in Eq. (25). 


Eff total = 


Power net + mu 2 • LHVu 2 
mf ee d • LHV feed 


(25) 


For all the case studies, we calculate the revenue in a simple 
way that only takes into account the raw material costs variation 
and the profit obtained by selling all of the electricity and hydro¬ 
gen produced. Therefore, for comparison purposes, we assume that 
investment in the plant and its operation and maintenance costs 
are the same for all case studies. The assumed prices for each item 
are: coal 45EUR/t, petcoke 75EUR/t, orujillo 65EUR/t, electricity 
0.15 EUR/kW h and 3 EUR/kg for hydrogen. We assume 7200 work¬ 
ing hours per year. 

3. Results 

Different products are obtained when topological changes are 
made in the modelled flowsheet, according to the fraction of syngas 
that enters the CC or the C0 2 removal process: 

• All the syngas is used in the CC. 

• Co-production of power and H 2 . 

• Production of H 2 only. 

• Production of H 2 only. The PSA purge gas is recycled to be used 
in the CC. 

• The H 2 stream that is separated after the Rectisol system is used 
in the CC. 

The superstructure has a set of splitters that allows the afore¬ 
mentioned processes to be carried out, as represented in Fig. 3. In 
this study, we consider the alternatives marked with dotted green 
lines, color figures are available in the online version. Several input 
variables, such as flowrates and flow ratios, remain constant for all 
the topological options in all the case studies that have the same 
feedstock. These values are summarised in Table 5. This section is 
divided into three main parts, based on the different co-gasification 
and co-production studies performed with the model. Firstly, we 
study the performance of the base case working in IGCC mode (50% 
coal, 50% petcoke on a mass basis) by carrying out a sensitivity anal¬ 
ysis (SA) that takes into account the variation in both oxygen and 
steam flows to the gasifier. Parameters such as gasifier efficiency or 
global efficiency are considered as a matter of comparison. This base 
case is described in Section 3.1 . Secondly, binary and ternary blends 
are analysed and compared in terms of H 2 production in the gasi¬ 
fier block, the efficiency results and the raw gas composition. Then, 
co-gasification with co-production case studies are evaluated (IGCC 
mode, combined with hydrogen production). See Section 3.2 for a 
summary of these results. Finally, Pareto curves are considered for 
the case studies that focus on a multi-objective appraisal accord¬ 
ing to environmental, economic and technical criteria. These are 
described in Section 3.3. The discussion of model trends in Sections 
3.1 and 3.2 are developed by performing various SAs. Differences 
between each case study lie in the raw material composition and 
changes in the oxygen ratios for the gasifier and the combustion 
chamber to adjust T gasi j and TOT (according to the technical data in 
Table 1). SAs are performed for each change in gasifier feedstock: 
the ER is changed in an interval of ±30% the base case value, while 
the steam ratio is varied between 0.15 and 0.4 for all case studies. 
Eight values are considered in both ranges. Accordingly, 64 points 
are given in each analysis. The CC-fed oxygen ratio is kept constant. 
In the case studies reported in Table 5, we consider the extreme 
cases of production of electricity or hydrogen for the market. 

3. I . Coal and petcoke co-gasification for electricity production 

To evaluate the influence of process variables on trends in model 
output variables, the results obtained from the SAs by varying the 
two gasifier ratios are represented in two-dimensional graphs, in 
which the axes correspond to the input variables and the output 
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Table 5 

Co-gasification of coal, petcoke and orujillo main results. 



Cl 

C2 

C3 

C4 

C5 

C6 

C7 

C8 

C9 

Inputs 

Coal (mass fraction) 

0.5 

1 

0 

0 

0.45 

0.35 

0.25 

0.15 

0.05 

Petcoke (mass fraction) 

0.5 

0 

1 

0 

0.45 

0.35 

0.25 

0.15 

0.05 

Orujillo (mass fraction) 

0 

0 

0 

1 

0.1 

0.3 

0.5 

0.7 

0.9 

ER 

0.422 

0.520 

0.370 

0.406 

0.420 

0.418 

0.414 

0.411 

0.408 

H 2 0 ratio (from gasifier) 

0.16 

0.16 

0.16 

0.16 

0.16 

0.16 

0.16 

0.16 

0.16 

Outputs 

LHV rg (MJ/kg) 

9.5 

6.6 

11.4 

7.9 

9.4 

9.1 

8.8 

8.4 

8.1 

CGE (%) 

75.2 

59.9 

82.9 

60.7 

74.0 

71.1 

68.3 

65.3 

62.3 

Liquid sulphur (t/h) 

3.32 

0.98 

5.74 

0.11 

3.00 

2.37 

1.73 

1.08 

0.43 

Electricity production (fraction of syngas to H 2 
0 2 ratio CC 2.45 

= 0) 

2.13 

2.54 

2.45 

2.48 

2.46 

2.46 

2.45 

2.44 

Effcc(%) 

57.7 

58.3 

57.5 

60.8 

57.9 

58.5 

59.1 

59.7 

60.4 

Net power (MW) 

285.3 

151.9 

423.2 

243.9 

281.6 

273.6 

265.3 

256.9 

248.4 

Effgiotoi (%) 

42.1 

34.3 

45.7 

36.8 

41.6 

40.6 

39.6 

38.5 

37.4 

CO 2 in flue gas (t/h) 

223.2 

131.8 

317.9 

170.4 

217.9 

207.2 

196.6 

186.1 

175.6 

Hydrogen production (fraction of syngas to H 2 
EffcCapture (%) 68.9 

= 1) 

69.4 

61.3 

70.7 

69.0 

69.3 

69.6 

69.9 

70.4 

H 2 (t/h) 

10.4 

5.5 

13.6 

8.6 

10.1 

9.7 

9.4 

9.0 

8.7 

Effglobal H2 (%) 

50.2 

40.8 

50.9 

42.8 

49.6 

48.1 

46.6 

45.1 

43.6 

Net power (MW) 

-91.9 

-43.1 

-120.9 

-35.2 

-75.3 

-52.1 

-47.3 

-42.2 

-43.5 

C0 2 to final disposal (t/h) 

207.3 

124.9 

277.6 

158.5 

202.6 

191.9 

180.8 

171.6 

155.1 


variable changes are plotted using contour lines. Blue points on the 
graphs are the results of SAs, while the red star point is the highest 
value, of the output variable in the SA. The base case SA considers 
an oxygen ratio variation from 0.3 to 0.55. 


Fig. 7 reports the variation observed for T gasi f, TOT and the LHV rg . 
The steam’s role as a moderator in the gasifier is clearly shown 
(more steam, leads to lower temperatures). In addition, more inlet 
oxygen is associated with higher temperatures. The SA shows gasi- 



Fig. 7. Iso-lines for T gasi f, TOT, and LHV rg in ° C and J/kg respectively, for the base case. 
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Mole h 2 



Fig. 8. Iso-lines representing the molar fraction (%) 


fication temperatures in the range of 800° C and 2400 °C, even 
when the gasification zone is in the middle and limited by CH 4 and 
sulphur oxide production. The TOT is between 550 0 C and 490 °C. In 
the modelled installation, temperatures below 540 0 C are not feasi¬ 
ble in operating regions of the HRSG, due to the steam temperature 
requirements. As regards the LHV rg , the highest values are obtained 


Mole CO 



ER 


of raw gas components in dry basis for the base case. 

for the combination of the lowest H 2 0 ratio with the lowest ER. Syn¬ 
gas compositions (H 2 , CO, C0 2 , and CH 4 ) are represented in Fig. 8. 
The occurrence of CH 4 is mainly driven by the 0 2 amount coupled 
with high gasifier temperatures (see Fig. 7, left and Fig. 8, bottom 
right). Non-linear shapes are found for H 2 , CO and C0 2 . The highest 
H 2 production is found for an ER close to 0.4 and for the highest 



Fig. 9. Iso-lines for gas and steam cycles power production (MW) for the base case. 










































































M. Perez-Fortes et al. / Computers and Chemical Engineering 35 (2011) 1501-1520 


1513 



Fig. 10. Iso-lines for CGE (Eq. 16), Effcc (Eq. 18) and Eff g i obal (Eq. 19) in % for the base case. 


ratio of steam that is fed into the gasifier, which is in clear agree¬ 
ment with the endothermic Eqs. (5) and (6). A clear relationship is 
found between C0 2 and H 2 composition, but the highest values of 
C0 2 are observed for the maximum amount of 0 2 and steam. As 
expected, the CO behaviour is the opposite of the C0 2 behaviour. 
Figs. 9 and 10 are highly interrelated. Firstly, ST and GT power is 
plotted. The amount of power produced by the GT reaches a max¬ 
imum when the ER and the steam ratio are at the lowest levels. In 
contrast, ST production is greater at higher oxygen ratios but lower 
steam ratios. In both cases, the 0 2 ratio has a greater influence than 
the steam ratio. The maximum GT power production coincides with 
the maximum level of CO production. In turn, the ST power produc¬ 
tion is larger when the consumption of steam in the system is at 
the lowest values. 

As shown in Fig. 10 (top left), the highest value of CGE (defined in 
Eq. (16)) coincides with that obtained for the maximum GT power 
production. This is in concordance with the higher LHV rg values 
found in Fig. 7, given that the fuel LHV remains constant. The Eff g i oba i 
shows the same behaviour due to the fact that most of the power 
comes from the GT and the fuel LHV remains constant. The highest 
value for Eff cc coincides with the highest ST power production and 
therefore with the highest gasification temperatures, at which the 
WHB recovers more heat for steam production. 

The ER has a greater influence on gasification than the steam 
ratio. CGE is the main parameter that influences the final global 
efficiency of the plant when it is working in IGCC mode, since it is 
the GT that produces more power. 


3.2. Coal, petcoke and biomass co-gasification for the 
co-production ofH 2 and electricity 

In this section, we consider the use of coal, petcoke and orujillo 
individually and in mixtures. Firstly, we assess the influence of the 
feedstock on gasification and on the IGCC mode of operation by 
changing the oxygen and steam ratios. Secondly, all combinations 
are studied by taking into account 100% H 2 or electricity produc¬ 
tion. Table 5 shows the different mixtures that are considered 
and the most significant model output values for both modes of 
working. Thirdly, several SAs are performed to evaluate the impact 
of feedstock mixtures and the split fraction of syngas to produce 
electricity or H 2 . 

In Fig. 11, we compare the H 2 production in the gasifier when 
it is fed by pure raw materials. In all cases, we observe a similar 
behaviour to that found in the base case, in which the maximum 
amount of H 2 is produced for the highest ratios of H 2 0 (see Fig. 8). 
Orujillo, which has the highest proportion of oxygen and hydro¬ 
gen in its composition, achieves the highest proportion of H 2 with 
the lowest ER values. It is followed by petcoke, which has a very 
similar value. Coal leads to a higher ER. Figs. 12 and 13 compare 
ternary blends in terms of CGE and Effcc • Firstly, the overall trend in 
CGE results is the same for all the blends, even though the absolute 
value is higher for mixtures with lower proportions of biomass. This 
is due to the effect of petcoke, which is the fuel component with the 
highest calorific value. In addition, the effect of the H 2 0 ratio is only 
noticeable at lower ER values, at which lower steam ratio values 
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100% Coal Mole H 2 100% Petcoke Mole H, 



0.3 0.35 0.4 0.45 0.5 

ER 

Fig. 11. Iso-lines for H 2 production (molar fraction % in dry basis) of the three pure components. 


imply higher CGEs. The tested ER interval variation for gasification 
is quite similar in all the graphs. Secondly, Effcc has slightly higher 
values when higher percentages of biomass are considered. For this 
efficiency value, the steam ratio does not have any significant influ¬ 
ence, as shown by the parallel isolines. Note that CGE and Effcc still 
follow opposite trends when ER is increased, presenting a trade-off 
that needs to be solved in order to set an ER value that maximizes 
Eff global- Higher values of Effcc are associated with higher values of 
ER (marked with red stars in Fig. 13). The addition of biomass to the 
base case, while the proportion of coal and petcoke is maintained, 
is shown in Fig. 14. A similar analysis was performed by Hernandez, 
Aranda-Almansa, and Serrano (2010), who evaluated the effect of 
biomass content on the syngas composition ratios H 2 /CO, CH 4 /H 2 , 
H 2 /C0 2 and C0/C0 2 , as well as on CO and H 2 content and the CGE. 
We have added the CC efficiency, the whole plant efficiency and 
other syngas components to these outputs. Our results are obtained 
at higher temperatures than those used in the aforementioned 
paper, and are obtained through an equilibrium Gibbs reactor. Both 
of these factors result in different trends. The main reason for the 
discrepancies may be the fact that biomass synergies are not appre¬ 
ciated in our model, and that the syngas compositions reported in 
the aforementioned paper might not reach chemical equilibrium. 
To sum up, the addition of biomass enhances the production of 
hydrogen, but leads to lower CGE which, in turn, brings about an 
overall reduction in efficiency. Regarding the syngas composition, 
the increase in biomass entails less CO and more C0 2 . 

Secondly, the reference used in all scenarios is the plant feed 
flowrate, or in other words, the maintenance of the load of the 


plant at approximately 110 tons/h. Hence, the final plant output is 
different in all case studies. The section outputs in Table 5 refers to 
the main parameters selected after gasification and after the syn¬ 
gas cleaning units, which remain the same for both final syngas 
applications. The LHV of the syngas just after the gasifier ( LHV rg ) 
is different for each feed. It is highest for the petcoke, followed by 
the base case, the orujillo and the coal scenario. For C6-C9, it varies 
in proportion with the amount of biomass introduced. It is directly 
related to the LHV of the raw materials, as occurs with the CGE. The 
liquid sulphur produced depends on the feedstock composition, 
and therefore on the sulphur percentage according to the ultimate 
analyses. The sulphur release efficiency is approximately 96% in all 
case studies. The final value of C0 2 emissions to the atmosphere 
(after the syngas combustion in the CC) is also represented. This is 
in agreement with the C composition of the feedstock. Therefore, 
the petcoke emits more C0 2 in absolute values. Hydrogen pro¬ 
duction is represented by its efficiencies and its final streams. The 
efficiencies follow an analogous trend to the CC option tendency: 
Eff g iobai H2 diminishes when the biomass fraction in the feedstock 
mixture increases. Again, the petcoke performs best, as it produces 
more H 2 with the highest efficiency. The final amount of hydrogen 
is not proportional to the final amount of hydrogen produced in 
the gasifier. This is due to the CO composition: the final amount 
of hydrogen recovered depends on the H 2 /CO ratio before the WGS 
reaction. Therefore, this ratio should be optimised in carbon capture 
applications, rather than the value of H 2 alone. 

The flowrate of the captured C0 2 (which is sent to geological 
storage) is also reported. This stream differs from the C0 2 stream 
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10% Biomass CGE 30% Biomass CGE 



50% Biomass CGE 



ER 


70% Biomass CGE 



90% Biomass CGE 


100% Biomass CGE 



Fig. 12. Iso-lines for CGE values (%) (Eq. (16)) of the different biomass blends. 


from the flue gas in the final amount of C0 2 that is emitted to the 
atmosphere through the purge PSA gas. Case studies with higher 
C0 2 emissions result in scenarios with lower ratios of carbon to car¬ 
bon storage in the flue gas. Again, this is due to the WGS reactors: 
the equilibrium constants of the two reactors at different temper¬ 
atures depend on the composition of the reactants. Therefore, the 
syngas composition is a critical parameter, as observed. 

To sum up, the feedstock compositions that produce the most 
electricity also produce the most hydrogen. In addition, the H 2 /CO 
ratio should be optimised in the final carbon capture application. 


Four SA with feedstock and product variations are performed. 
3D graphs are drawn up here to represent the variation in the cho¬ 
sen variables, in terms of the split fraction of syngas that goes to H 2 
production and the fraction of biomass (orujillo) introduced into 
the 50:50 coal and petcoke mixture. The representative variables 
that best describe the co-production scenarios (see Fig. 15) are the 
net power produced when the plant is working in IGCC mode, the 
amount of MW of H 2 produced when it is working with carbon 
removal units, and the Eff tota i (Eq. 25), which takes into account the 
efficiency of the plant and evaluates both, the products and the rev- 
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10% Biomass Eff^ 30% Biomass Eff^ 



Fig. 13. Iso-lines for the CC efficiency evaluation (%) (Eq. (18)) of the different biomass blends. 


enue, as explained in Section 2.3. The top-left graph of Fig. 15 shows 
that as the fraction of syngas to hydrogen increases, the power 
obtained decreases until negative values are reached. At these val¬ 
ues, all the syngas produces H 2 and, as a result, the power produced 
by the ST only (which works using steam from the WHB) is not 
enough for the self-supply of electricity. The analysis shows that 
higher fractions of biomass, produce less net power. The highest net 
power is obtained when both split fractions are zero. The top-right 
graph shows the opposite behaviour, in feedstock terms. In this 


case, the highest H 2 energy flowrate is obtained when the biomass 
fraction is 0 and the H 2 split fraction is 1. If we match both perfor¬ 
mances, the bottom-left graph displays a behaviour that is more or 
less stable compared to the previous graphs: Eff tota i has a value that 
ranges between 35% and 43%. The highest values are found for the 
lowest fractions of biomass in the feedstock mixture, and for the 
highest fractions of H 2 produced. Finally, for the assumed prices, 
the bottom-right graph seems to show the opposite behaviour to 
that observed in the previous graph: even when the efficiency is not 
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Fig. 14. Effect of biomass addition in raw gas composition and in efficiencies. 


at its maximum, scenarios with higher production of electricity are 
favoured. 

As can be deduced from Fig. 16, a price-based indicator can be 
very subjective, since changes in prices can alter completely the 


selection of which case study is best. For instance, in these graphs 
we can observe how the revenue plot changes if the prices are modi¬ 
fied. The top-left graph shows the behaviour of the scenarios when 
the price of the electricity is reduced to 0.05EUR/kWh. The plot 
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Fig. 15. Analysis of the effect of biomass fraction variation with co-production of electricity and hydrogen. 
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Fig. 16. Effect of the price on the revenue behavior. Top-left graph has a decrease on the electricity price. Top-right graph, on H 2 price. Down-left and down-right graphs 
vary the biomass price: first, it is set to zero; then, we assume that we are paid for collecting it. 


follows the same trend as the Eff tota i plot and favours the hydro¬ 
gen production case studies. The top-right graph uses an H 2 price 
of 1 EUR/kg. In this case, the best behaviour is shown by the same 
scenarios as in the bottom-right graph in Fig. 15, favouring the use 
of no biomass and with the production of electricity alone. When 
the biomass price is zero, as in the bottom-left graph, electricity 
production is favoured at any biomass fraction. In the bottom- 
right graph, we are paid 75 EUR/t for collecting the biomass. In 
this case, the optimal scenario is the one with the highest biomass 
fraction and the production of 100% electricity. To sum up, the 
choice of an appropriate mixture of feedstocks and the amount 
of products is dependent on many criteria and is subject to the 
decision-maker’s interests. Nevertheless, this co-gasification and 
co-production analysis highlights the fact that the production of 
H 2 together with electricity is a more efficient process, in terms of 
Efftotai- However, the user has to consider that the Ejf g i oba i is reduced 
when carbon capture units are added to an already existing IGCC 
power plant. Accordingly, the feedstock should be increased to pro¬ 
duce the same amount of electricity as before. Again, this situation 
shows the importance of each case study and each specific context, 
such as the demands to be met, or the raw material and product 
prices. Hence, systematic tools are needed for the decision-making 
process, such as that developed in Section 3.3. 

3.3. Overall approach: Pareto curves 

The generation of Pareto frontiers is one of the many ways to 
aid in multicriteria decision-making problems. The evaluation cri¬ 
teria for different case studies include a set of model outputs that 
take into account a mixture of co-production parameters, thereby 
considering a multi-objective approach which includes the gasifier 
efficiency (CGE), the raw material cost, the C0 2 relative emissions in 
IGCC mode, the net power produced in IGCC mode and the hydro¬ 


gen (in MW) produced when carbon capture units are used. Fig. 17 
shows the different solutions in two-dimensional space. The Pareto 
frontier in each graph is the collection of case studies results that 
are not dominated by other configurations. Given that five indi¬ 
cators are considered in this case, then ten Pareto curves can be 
obtained by binary criteria comparison. 

In our case, we want to obtain the most appropriate feedstock 
composition to maximise the production of H 2 to be sold to the 
market as pure H 2 or to produce power. Therefore, the objective for 
each indicator is as follows: low raw material cost, low C0 2 relative 
emissions, high CGE and a high value of H 2 (in MW) or high pro¬ 
duction of net power. These parameters are used in the calculation, 
with the same working hours and prices as in Section 2.3. 

An analysis of the simulation results for each case study shows 
that C3 is the best scenario for maximising hydrogen production, 
net power and CGE. There is a trade-off between the raw material 
costs and the relative C0 2 emissions and all the other parameters. 
Consequently, only seven out of the ten possible Pareto curves are 
represented. 

For each Pareto curve in Fig. 17: 

• If CGE and raw material costs are considered, then the frontier is 
determined by C3-C1-C2. These are the base case and the case 
studies that use only coal and petcoke as raw materials. 

• In the case of C0 2 emissions with raw material costs the Pareto 
frontier contains all the scenarios except C3. Consequently, for 
the same (or lower) amount of relative carbon emissions there 
are cheaper options than petcoke. 

• When raw material costs are compared with the hydrogen and 
power produced, the result is the same Pareto frontier made up 
of the C3-C1-C2 scenarios. 

• When C0 2 emissions are compared with the net power and 
with the CGE, the result is a Pareto frontier formed by 
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Fig. 17. Pareto curves that evaluate changes in raw material composition. 


C4-C9-C8-C7-C3. The highest producers of power are not nec¬ 
essarily the highest producers of emissions related to the GWh 
produced. Moreover, as already stated, optimal CGE leads to opti¬ 
mal overall results in terms of energy produced. 

• C0 2 emissions vs. the hydrogen produced, results in a Pareto fron¬ 
tier made up of C4-C9-C8-C7-C3; the same frontier as that found 
in the previous case. 

Overall, C3 is always on the Pareto frontiers, except in the com¬ 
parison of carbon emissions with raw material cost. In contrast, Cl 


and C2 are favoured by the raw material cost. Although they have 
low values of CGE, net power produced and hydrogen produced, 
they are the cheapest options. Scenarios with 100% of biomass and 
the highest proportions of it (C7-C8-C9) are favoured according 
to the criteria defined in these Pareto frontiers. Therefore, for the 
purposes outlined here, the best options involve co-gasification 
with biomass, which is either used alone or with petcoke as the 
raw material. Coal and co-gasification of coal and petcoke are good 
choices regarding as they are economical. As is clearly pointed out 
here, the financial aspects follow market laws, which are not nec- 
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essarily governed by the most efficient choices. Moreover, changes 
in prices will change the Pareto frontiers. 

4. Conclusions 

The results provided by this validated IGCC model show its pre¬ 
diction capacity. The procedure involves evaluation scenarios that 
ease the decision-making task and take into account multiple objec¬ 
tives. The model can be used for changes in raw materials and 
in operation units. Comparisons indicate that the model’s results 
agree to a high degree with real power plant data. 

Oxygen has been shown to have a greater influence on gasifica¬ 
tion than steam. Biomass gasification and biomass co-gasification 
with coal and petcoke are attractive choices for climate change mit¬ 
igation and traditional fuels independence. The simulation results 
analysed by means of Pareto curves show that biomass is a suitable 
choice for net power and H 2 production, as it offers higher efficien¬ 
cies and lower carbon-related emissions. High CGE values lead to 
optimal overall results in terms of the electricity produced. Regard¬ 
ing prices of raw materials, its behavior associated to market laws 
does not reflect on the most efficient scenarios, which were not 
found as the most profitable ones. Concerning feedstocks and co¬ 
production, the feedstock compositions that lead to the most elec¬ 
tricity production also offer the best hydrogen production results. 

The objective of this approach was to conceive a model and a 
useful tool for finding optimal solutions that enable plant condi¬ 
tions to be improved. Depending on the final interests, this tool can 
help in the formulation of design options for planned changes or be 
used to assess trends in specific parameters. 
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